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We present a high-quality description of the deuteron electromagnetic form factors in a soft-wall 
AdS/QCD approach. We first propose an effective action describing the dynamics of the deuteron 
in the presence of an external vector field. Based on this action the deuteron electromagnetic form 
factors are calculated, displaying the correct 1/Q 10 power scaling for large Q 2 values. This finding 
is consistent with quark counting rules and the earlier observation that this result holds in confining 
gauge/gravity duals. The Q 2 dependence of the deuteron form factors is defined by a single and 
universal scale parameter k, which is fixed from data. 
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The experimental and theoretical study of the deuteron is one of the main focuses of hadronic physics during the last 
decades (for detailed reviews see e.g. Refs. 0 - 0 )- Many theoretical approaches have been applied to the problem of 
the deuteron form factors: per turbative QCD, chiral effective and phenomenological approaches, potential and quark 
models (see e.g. Refs. (lj|-|22|). For example, in potential models the nonrelativistic impuls approximation |3, [§] 
was used. It leads to deuteron form factors factorized in terms of the isoscalar combinations of the nucleon form 
factors. These approaches are able to describe data up to 0.5 GeV 2 , but deviate from data for higher Q 2 and are 
not consistent with quark counting rules. To include relativistic effects different types of relativistic nuclear models 
have been developed. One possibility is based on taking into account relativistic corrections in a u/c expansion of the 
nonrelativistic current (leading to so-called two-body interaction current diagrams) Such approaches are limited 

in their validity of the description of data up to 1-2 GeV 2 . There is a group of models based on relativistic Hamiltonian 
constraint dynamics, which uses certain phenomenological potentials (Argonne, Nijmegen, etc.) and three forms of 
quantization procedures (point, instant or front form) [see e.g. Refs. 12], Ll_3j] • Field-theoretical methods formulated 
in terms of hadronic (mesons, nucleons, A-isobars) degrees of freedom are used in a wide range of approaches. These 
include models based on the solution of a quasipotential 141 or on Bethe-Salpeter |7(| equations. These methods 
also include field theories quantized on the light cone 17. 18j. phenomenological Lagrangian approaches Hi and 
effective field theories treating the long-range dynamics explicitly while parametrizing the short-distance effects by 
contact interactions [for recent applications to deuteron form factors see e.g. Refs. [loj]]. Another class of approaches 
supposes to treat the deuteron in terms of fundamental degrees of freedom - quarks and gluons: nonrelativistic quark 
models 21, 2f| and perturbative QCD [g]. The analysis of Ref. [@j results in a prediction for the asymptotic large- 
momentum-transfer behavior of the deuteron form factors and the form of the deuteron distribution amplitude at 
short distances. Later on in Ref. [23| it was shown that field theories based on gauge/gravity duality, as proposed in 
Refs. |2J], produce the correct power scaling of hadronic form factors at large momentum transfer. This finding is 
consistent with the quark counting rules. 

The main advantage of our approach is that it gives a description of the deuteron electromagnetic (EM) form 
factors in terms of a single dimensional parameter k with the correct power scaling 1 /Q 10 at large Q 2 as predicted by 
perturbative QCD. Our approach is constructed as a holographic dual to perturbative QCD. It gives a good starting 
point for studying more complicated many-body nuclear systems. 

Encouraged by this property we focus on soft-wall anti-de Sitter/quantum chromodynamics (AdS/QCD) 

It is a version of a bottom-up approach based on the correspondence of string theory in AdS space and conformal field 
theory (CFT) in physical space-time. The formalism in soft-wall AdS/QCD is based on an effective action involving 
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five-dimensional fields propagating in AdS space, which are dual to the deuteron and the electromagnetic field. We 
apply our formalism to the calculation of the EM form factors of the deuteron. The deuteron itself is simply considered 
as a proton-neutron bound state. 


Note that there are already some applications of AdS/CFT and AdS/QCD to different 
baryon matter at finite temperature and baryon number density [ 2 £j , cold nuclear matter 


potential [31] ; p meson condensation at finite isospin chemical 
atomic nuclei |34|, nuclear matter to strange matter transition 3j 
for baryon many-body systems [37] (for reviews see e.g. Refs. [31: 


problems in nuclear physics: 
[30(, baryon-char ge c hemical 
heavy 


'Otential [32|, holographic nuclear matter c _^ 

self-bound dense objects Q, and mean-field theory 

SI). 


Our approach is based on an effective action, which in terms of the AdS fields d M (x, z) and V M (x, z ), is dual to the 
Fock component contributing to the deuteron with twist r = 6 , and the electromagnetic field, respectively, is given 

by 


S = / d 4 xdze~^ z) 


z)F mn (x, z) - D m d? N (x, z)D M d N (x, z) - ic 2 F MN (x, z)d\ I {x, z)d N {x, z) 


+ TKI 2 e2A(2) dMpNK ( x ’ z ) (i D Kd ] M {x, z)d N (x, z) - d\f(x, z)iD K d N (x , z) + H.c.^ 
+ ^ M (x,z) (jjl 2 + f/(z)) d M (x,z) , 


( 1 ) 


where A(z) = log (R/z), F MN (x,z) = d M V N (x, z) — d N V M (x,z) is the stress tensor of the vector field V M (x,z), 
D m = d AI — ieV M (x, z) is the covariant derivative, n 2 R 2 = (A — 1)(A — 3) is the five-dimensional mass; R is the AdS 
radius, <p{z) = n 2 z 2 is the background dilaton field; A = r + L is the dimension of the d AI (x, z ) field; L is the orbital 
angular momentum, and Md is the deuteron mass. U(z) is the confinement potential with 


U{z) 


R 2 


Un 


( 2 ) 


where the constant Uq is fixed by the value of the deuteron mass. In the following we work in the axial gauge for 
both vector fields d z (x,z) = 0 and V z (x, z) = 0. In our consideration we have two free parameters: k and Uq (the 
latter only relevant for the description of the deuteron mass). As it will be shown later, the parameters c 2 and C 3 are 
constrained by normalization of the deuteron electromagnetic form factors. 

First we perform a Kaluza-Klein (KK) decomposition for the vector AdS held dual to the deuteron 


d ll (x , z) = exp 


■ if(z) - A(z) 
2 


J2 d n( X )®n{z) , 


( 3 ) 


where (x) is the tower of the KK fields dual to the deuteron fields with radial quantum number n and twist-dimension 
r = 6 , and & n (z) are their bulk profiles. 

Then we derive the Schrodinger-type equation of motion (EOM) for the bulk profile & n (z) with 


dz 2 


4(L + 4 ) 2 - 1 
4^2 


4 2 

■ AC Z 


k 2 U 0 


4?n(-Z) — 


( 4 ) 


The analytical solutions of this EOM read 


$n(z) 

Ml n 


2n\ 


(n + L + 4)! 


C L +5 _L+9/2 -k 2 2 2 /2 


L L n + \* 2 z 2 ), 


Ah 2 


n - 


L + 5 Uq 


( 5 ) 


where L™(x) are the generalized Laguerre polynomials. Restricting to the ground state (n = 0, L = 0) we get 
Md = 2 k \j\ + Using the central value of data for the deuteron mass Md = 1.875613 GeV and k = 190 MeV 
(fitted from data on electromagnetic deuteron form factors), we fix Uo = 87.4494. We can compare this value for the 
deuteron scale parameter to the analogous one of kn defining the nucleon properties - mass and electromagnetic form 
factors. In Ref. 0] we fixed the value to kn ~ 380 MeV, which is 2 times bigger than the deuteron scale parameter n. 
The difference between the nucleon and deuteron scale parameters can be related to the change of size of the hadronic 
systems - the deuteron as a two-nucleon bound state is 2 times larger than the nucleon. 
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In the case of the vector field dual to the electromagnetic field we perform a Fourier transform with respect to the 
Minkowski coordinate 


V»(x,z) = j ^e^VMV^z) 


where V ( q , z) is its bulk profile obeying the following EOM: 


d z 


o-v(z) 


o~<P(.z) 


d z V(q,z ) + g 2 - V(q,z) = 0. 


Its analytical solution [25] can be written in the form of an integral representation introduced in Ref. [42 

1 


V(Q,z) = k 2 z 2 J 


dx 


{l-xf 


0 -k, 2 z 2 x/( 1-x) a 


Q 


x ° > a = i^2’ Q =-Q ■ 


( 6 ) 


(7) 


( 8 ) 


The gauge-invariant matrix element describing the interaction of the deuteron with the external vector field (dual to 
the electromagnetic field) reads 


M ? nv (P’P') = -(^Gi(Q 2 )e*{p')-e(p) - C ^ 2 1 e*(p') ■ qe(p) ■ q'j (p + p'Y 
- G 2 (Q 2 ) (e<*(p) e*(p 1 ) • q - ?<*(]/) e(p) • g) 


(9) 


where e(e*) and p(p') are the polarization and four-momentum of the initial (final) deuteron, and q = p' — p is the 
momentum transfer. The three EM form factors Gi, 2,3 of the deuteron are related to the charge Gc , quadrupole Gq 
and magnetic Gm form factors by 

Gc = Gi + —TdGQ , Gm = G 2 , Gq = Gi — G 2 + (1 + Td)Gs, t d = . ( 10 ) 

These form factors are normalized at zero recoil as 

Gc(0) = 1, Gq( 0 ) = M 2 Qd = 25.83, G M (0) = = 1.714, ( 11 ) 

Mn 

where Md and Mn are deuteron and nucleon masses, and Qd = 7.3424 GeV " 2 and p d = 0.8574 are the quadrupole 
and magnetic moments of the deuteron. Since the deuteron is a spin-1 particle it has three EM form factors in the 
one-photon-exchange approximation, due to current conservation and the P and G invariance of the EM interaction. 

We illustrate the algorithm for calculating the deuteron form factors, considering a particular case of the form factor 
Gi(Q 2 ), which is generated by the second term in the effective action <[lj), 

= J d 4 xdz e~^ z ) eVfj,(x, z) z)d l '(x, z) — dj,(x, z)id^d v (x, z)j . (12) 

Next we use the Kaluza-Klein decomposition (0 for the five-dimensional fields d v {x,z) and dl(x,z) (restricting to 
the contribution of the ground states with n = 0) and perform the Fourier transform for d„(x), dj,(x), 

dAx) = / WY e ~ ipx eAp) > 4 (*) = / (13) 

and V/j.(x,z) [see Eq. ©]■ Substituting expressions © and m in action m and integrating over x and z, we get 



where 1 (p, p') is part of the invariant matrix element of the d + 7 —> d transition containing the contribution of 

the form factor Gi(Q 2 ) 

MP^\p,p') = -{p+p'Y e*{p') ■e{p)G 1 {Q 2 ). ( 15 ) 
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In our approach the deuteron form factor G\(Q 2 ) = F(Q 2 ), where F(Q 2 ) is the twist -6 hadronic form factor, which 
is given by the overlap of the square of the bulk profile dual to the deuteron wave function (twist -6 hadronic wave 
function) and the confined electromagnetic current 

OO 

F(Q 2 ) = J dz <t> 2 (z) V(Q, z) = ( 16 ) 

o 

where a = Q 2 /{An 2 ). This formula follows from the general and universal formula for the hadronic form factor with 
twist r derived in Ref. [27| in terms of the bulk profile <j>A z ) = yj (r- 2 )! z T ~ 3 l 2 e~ K 2 / 2 dual to the hadronic wave 
function with twist r: 

OO 

FAQ 2 ) = j dz <f>A z ) V(Q, z) = ( 17 ) 

0 

Therefore, Eq. (fTBl) is the particular case of Eq. 1171) for r = 6 . 

By analogy we calculate the other two deuteron form factors Gi and G 3 , which are expressed in terms of the same 
universal factor F(Q 2 ): 


Gi(Q 2 ) = CiF(Q 2 ), * = 2,3. 

The parameters C 2 and C 3 are defined by normalization of the deuteron form factors as: 

c 2 = G m ( 0) = 1.714, c 3 = G m (0) + G q ( 0 ) - 1 = 26.544. 


(18) 


(19) 


Note that the form factor F(Q 2 ) has the correct power-scaling F{Q 2 ) ~ 1/(Q 2 ) 5 at large Q 2 —> 00 . It can also be 
written in the Brodsky-Ji-Lepage form derived within perturbative QCD. The deuteron form factor is factorized in 
terms of the nucleon form factor Fn(Q 2 /A) and the so-called “reduced” nuclear form factor fd{Q 2 ) 0: Fd(Q 2 ) = 
fd(Q 2 ) F^(Q 2 / 4). Our result reads 


F d (Q 2 ) = F(Q 2 ) 


T(6)T(a + l) 

T(a + 6) 


_5!_ 

(a + 1)... (a + 5) 


fd(Q 2 )F 2 (Q 2 / 4) 


( 20 ) 


where our predictions for fd{Q 2 ) and Fn(Q 2 /A) are 


f d (Q 2 ) 


30(q + l)(o + 2) 

(a + 3)(a + 4)(q + 5) 


Fn(Q 2 / 4) 


2 

(q + l)(q + 2 ) 


( 21 ) 


where a = Q 2 /(Ak 2 ). Our results for the charge Gc{Q 2 ), quadrupole Gq(Q 2 ) and magnetic Gm(Q 2 ) form factors are 
shown in Figs.1-3. The shaded band corresponds to values of the scale parameter k in the range of 150 MeV < k < 
250 MeV. An increase of the parameter k leads to an enhancement of the form factors. The best description of the 
data on the deuteron form factors is obtained for k = 190 MeV and is shown by the solid line. Data points are taken 
from Refs. To quantify the quality of the fit with k = 190 MeV we indicate the y 2 values for the three deuteron 

form factors: y 2 = 0.2 for (Gc), X 2 = 13-8 for (Gq) and x 2 = 2.3 for (Gm )■ We would also like to point out that with 

k = 190 MeV our result for the deuteron charge radius rc = (—6dGc{Q 2 )/dQ 2 ^^) 1 / 2 = yj — Qd = 1-846 fm 
compares well with data, rc = 2.130 ± 0.010 fm [l|. 

In conclusion we stress again the main result of this paper. Using the soft-wall AdS/QCD model we calculate the 
deuteron electromagnetic form factors, which are given by analytical expressions in terms of a universal twist -6 form 
factor F(Q 2 ) relevant for the deuteron — hadronic system with six partons. Our framework gives a description of the 
deuteron in terms of two free parameters — the dimensional parameter k and the confinement parameter Uq. The 
parameter k is fixed by the scale of the deuteron form factors and the parameter Go is fixed through the deuteron 
mass. 
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FIG. 1: Charge deuteron form factor. 





FIG. 2: Quadrupole deuteron form factor. 





FIG. 3: Magnetic deuteron form factor. 










